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Abstract 

We consider a supersymmetric version of the standard model extended by an additional U{1)b-l- 
This model can be embedded in an mSUGRA-inspired model where the mass parameters of the 
scalars and gauginos unify at the scale of grand unification. In this class of models the renormaliza- 
tion group equation evolution of gauge couplings as well as of the soft SUSY-breaking parameters 
require the proper treatment of gauge kinetic mixing. We first show that this has a profound 
impact on the phenomenolgy of the Z' and as a consequence the current LHC bounds on its mass 
are reduced significantly from about 1970 GeV to 1790 GeV. They are even further reduced if 
the Z' can decay into supersymmetric particles. Secondly, we show that in this way sleptons can 
be produced at the LHC in the 14 TeV phase with masses of several hundred GeV. In the case 
of squark and gluino masses in the multi-TeV range, this might become an important discovery 
channel for sleptons up to 800 GeV (900 GeV) for an integrated luminosity of 100 fb~^ (300 fb~^). 
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I. INTRODUCTION 



The LHC is rapidly extending our knowledge of the TeV scale. As there are currently no 
signs of new physics, models beyond the standard model (SM) have begun to be severely 
constrained. One of the most popular model classes is that of supersymmetric extensions, 
in particular the minimal supersymmetric standard model (MSSM). As the MSSM itself 
has over 100 free parameters, mainly models with a smaller set of parameters, like the 
constrained MSSM (CMSSM) with five parameters [HE], are studied. This model has been 
fitted to various combinations of experimental data: see e.g. |3HTT]. indicating that gluinos 
and squarks have masses in the multi-TeV range. In particular, if one wants to explain the 
recent hints of a SM-like Higgs boson with mass of 125 GeV fn[ [TB] , the CMSSM becomes 
more and more unlikely pLlj- However, this is mainly due to the strong correlations between 
the various masses of the supersymmetric particles. Once these are given up, a light SUSY 
spectrum is still compatible with LHC data as discussed in [TH - IT8] . 

In the framework of constrained models, a way out of this tension is to consider extended 
models such as the next-to-minimal supersymmetric standard model (NMSSM) [191 120] . 
Indeed, it has been shown that even in the constrained NMSSM a Higgs mass of 125 GeV 
can be explained [2T] 122] and in the generalized NMSSM these masses are obtained with 
even less fine-tuning [23] • A second possibility is to consider models with extended gauge 
structures, since then the upper bound on the lightest Higgs boson is also relaxed due to 
additional tree-level contributions [2H - I29] . Such models arise naturally in the context of 
embedding the SM gauge group in a larger group such as 5*0(10) or Eq: see e.g. [30ti35] . 
which can also nicely explain neutrino data via the seesaw mechanism [36ti39] . It can be 
shown that such models can have a Z' with a mass in the TeV range [lOj HI] . Z' searches 
are therefore among the targets of the Tevatron and LHC collaborations and bounds on its 
mass have been set [12HI5] . These bounds depend on the concrete gauge group and the 
couplings of the Z' to the SM fermions: see e.g. |46l450] . For reviews on various Z' models 
see e.g. [511152]. 

It has been known for some time that decays into supersymmetric particles can also 
strongly impact the phenomenology of the Z' [531 . The LHC phenomenology of super- 
symmetric U{1) extensions has been discussed in the context of S'O(IO) and Eq embeddings 
in [551457] and in a more general class of models with a U{1)b-xL in [HE]. In Ref. [56] a 
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U{1)b^l extension of the MSSM is also discussed. Furthermore, it has been pointed out 
that sleptons, charginos, and neutrahnos can be directly produced via a heavy Z' at the 
LHC even if they have masses of hundreds of GeV. 

f/(l) extensions of the SM have a peculiar feature, namely, gauge kinetic mixing [59V - 
1^ . In the previous studies either it has been argued that these effects are small or they 
have been completely ignored. However, it has also been pointed out in [521 [S3] that gauge 
kinetic mixing can be important in Z' searches in the case of Eq embeddings. Moreover, 
it has recently been shown that gauge kinetic mixing is important for the spectrum in 
U{1)b-l extensions of the MSSM [6l|. Such models can be embedded, for example, in a 
string-inspired Eg x Eq gauge group [651467] . 

In this paper we will explicitly show how existing collider bounds on the Z' mass are 
changed once gauge kinetic mixing is taken into account. We find that this can reduce the 
bounds by about 200 GeV, as the couplings of Z' to the SM fermions depend on this mixing. 
This clearly affects also the cross section for SUSY particles produced via a Z'. However, 
as we will demonstrate, LHC with ^/s = 14 TeV and a luminosity of 100 fb~^ should be 
able to discover sleptons via a Z' with slepton masses up to 800 GeV, provided the Z' is not 
much heavier than 2.8 TeV. In the case of 300 fb~^, this extends to slepton masses of about 
900 GeV up to Z' masses of 3.1 TeV. 

The remainder of this paper is organized as follows: In Sec. |Tl] we briefly summarize 
the main features of the model and its particle content. In Sec. |III| we first discuss how 
gauge kinetic mixing as well as supersymmetric final states affect the Z' phenomenology in 
the context of constrained models. Afterwards we discuss the possibilities of the LHC to 
discover SUSY particles, in particular charged sleptons, via a Z' . Here we will depart from 
the universality assumption as this mainly depends on the slepton and Z' masses. Finally 
we draw our conclusions in Sec. |IV[ In the Appendix we give the couplings of the Z' to the 
scalars and fermions, including terms arising due to gauge kinetic mixing effects. 

II. THE MODEL 

In this section we discuss briefly the particle content and the superpotential of the model 
under consideration and we give the tree-level masses and mixings of the particles important 
to our studies. For a detailed discussion of the masses of all particles as well as of the 
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TABLE I. Chiral superfields and their quantum numbers with respect to SU{'i)c ® SU{2)l ® 
U{1)y® U{1)b-l. 

corresponding one-loop corrections, we refer to In addition, we recall the main aspects 
of f/(l) kinetic mixing since it has important consequences for the phenomenology of the Z' 
within this model. 



A. Particle content and superpotential 

The model consists of three generations of matter particles including right-handed neu- 
trinos which can, for example, be embedded in 50(10) 16-plets [651467] . Moreover, below 
the grand unified theory (GUT) scale the usual MSSM Higgs doublets are present, as well 
as two fields 7] and f] responsible for the breaking of the U{1)b-l- The vacuum expectation 
value of 7] induces a Majorana mass term for the right-handed neutrinos. Thus we interpret 
the B — L charge of this field as its lepton number, and likewise for f/, and call these fields 
bileptons since they carry twice the usual lepton number. A summary of the quantum num- 
bers of the chiral superfields with respect to SU{?>)c x SU{2)l x U{1)y x U{X)b~l is given 
in Table [H 
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The superpotential is given by 

+ Y^v'rL,Hu -f/vfi +Y:^9::r^u^, (1) 
and we have the additional soft SUSY-breaking terms: 

^SB =^MSSM - Xb^b'^BB' - ]^\b'^B'^B' " "^^l??!^ " - '^I'^/iji^iT^'j 

- vvB,, + T^H^u^L, + T!^r^D^u] , (2) 

where i,j are generation indices. Without loss of generality one can take and -B^/ to be 
real. The extended gauge group breaks to SU{3)c^U{l)em as the Higgs fields and bileptons 
receive vacuum expectation values (vevs): 

^d=^ {(^d + Vd + i(j)d) , = ((T„ + w„ + (3) 

r/ = (d^ + t;^ + i(f)^) , f]= ^{afi + Vf^ + . (4) 

We define tan/3' = {v^/vfi) in analogy to tan/3 = {vu/vd) in the MSSM. 



B. Gauge kinetic mixing 

As already mentioned in the Introduction, the presence of two Abelian gauge groups in 
combination with the given particle content gives rise to a new effect absent in the MSSM 
or other SUSY models with just one Abelian gauge group: gauge kinetic mixing. In the 
Lagrangian, the combination 

- XaftF-^'^^F;,, a^h (5) 

of the field-strength tensors is allowed by gauge and Lorentz invariance [59] because F""'^" 
and F^'^^^ are gauge invariant quantities by themselves, 

Even if such a term is absent at a given scale, it will be induced by renormalization group 
equation (RGE) running. This can be seen most easily by inspecting the matrix of the 
anomalous dimension, which at one loop is given by 

where the indices a and b run over all U{1) groups and the trace runs over all fields with 
charge Q under the corresponding U{1) group. 



For our model we obtain in the basis (f/(l)y, U{1)b-l) 



7 = I ^ I AT, (7) 



167r2 
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and we see that there are sizable off-diagonal elements. N contains the GUT normalization 
of the two Abelian gauge groups. We will take as in ref. |6l] for U{1)y and for 
U{1)b-l, i-e. N = diag(-y|, -^1) and obtain 



The largeness of the off-diagonal terms indicate that sizable U{1) kinetic mixing terms are 
induced via RGE evaluation at lower scales, even if at the GUT scale they are zero. In 
practice it turns out that it is easier to work with noncanonical covariant derivatives instead 
of off-diagonal field-strength tensors such as in Eq. ([s]). However, both approaches are 
equivalent [HH]- Hence in the following we consider covariant derivatives of the form 

D, = d,-iQlGA, (9) 

where is a vector containing the charges of the field with respect to the two Abelian 
gauge groups, G is the gauge coupling matrix 

G = I I , (10) 

9BY QBE 

and A contains the gauge bosons A = {A^, A^)'^. 

As long as the two Abelian gauge groups are unbroken, the following change of basis is 
always possible: 

A = I ] ^ A'= i ] = RA, (11) 

where R is an orthogonal 2x2 matrix. This freedom can be used to choose a basis such that 
electroweak precision data can be accommodated easily. A particular convenient choice is 
the basis where Qby = because then only the Higgs doublets contribute to the entries in 
the gauge boson mass matrix of the U{1)y^ SU (2) ^ sector and the impact of rj and r/ is only 



in the off-diagonal elements as discussed in Sec. II D Therefore, we choose the following 
basis at the electroweak scale 



/ _gYYgBB — QybQby _ 
9yy / o — - 9i 



V fi'ls + 9%y 



9'bb =\/91b + 91y = 9bl 



9yb 



9yb9bb + 9by9yy 



91y 



9'by =0 • 



(12) 

(13) 
(14) 

(15) 



C. Tadpole equations 

Having in mind mSUGRA-like boundary conditions for the soft SUSY-breaking parame- 
ters, we solve the tadpole equations arising from the minimization conditions of the vacuum 
with respect to |yUp, 5^, and 5^/. Using = v"^ + f ? and f ^ = f J -|- f^, we obtain at 
tree level 

I/WP =\{{^99bl^'' cos(2/3') - Aml^ + 4m^J sec(2/3) - 4(m|^ + ml^j - [gl + f + gi^v"^^ 

(16) 

5^ = - ^ ( - '2g9BLX^ cos(2/3') + Aml^ - Am%^ + (^gf + f + gj^v^ cos(2/3)) tan(2/3) , 

(17) 

I/if =l{-2 (glLX^ + ml + mj) + {2ml " 2^ + 99blV^ cos(2/3)) sec(2/3')) , (18) 
B^, ( - 2glj^x^ cos(2/3') + 2m^ - 2m? + gg^^v^ cos(2/3)) tan(2/3') . (19) 

In the numerical evaluation we take also the one-loop corrections into account as discussed 
in [6l]. The phases of /x and /i' are not fixed via the tadpole equations and thus are taken 
as additional input parameters. As the phases are not important for our considerations, we 
set them to zero, e.g. sign(/i), sign(yu') > 0. 



D. Gauge boson mixing 



Because of the presence of the kinetic mixing terms, the B' boson mixes at tree level 
with the B and bosons. Requiring the conditions of Eqs. (12)-(15) means that the 



corresponding mass matrix reads, in the basis {B, W^, B') 



-Ir 



(20) 



In the hmit (7 — ?■ both sectors decouple and the upper 2x2 block is just the standard mass 
matrix of the neutral gauge bosons in electroweak symmetry breaking. This mass matrix 
can be diagonalized by a unitary mixing matrix to get the physical mass eigenstates 7, Z, 
and Z'. Because of the special form of this matrix, the corresponding rotation matrix can 
be expressed by two mixing angles Qw and 0'^^ as 



W 



v 



cos Qw — cos Q'w sin Qw sin Qw sin 9V 
sin Qw cos 9^^ cos Q'w — cos Qw sin Q'w 
sin V cos 6 V 



where Q'w can be approximated by [TO] 



tan 26' 



w 



1-g^9\^9l 



+ 16{x/vyg%^ ~ 9l- 9\ 



The exact eigenvalues of Eq. ( 20 ) are given by 



(21) 



(22) 



= 0, 



' = \ {{9\^ 9l^fV ^ ^9\l^''^ 



Expanding these formulas in powers of v /x\ we find up to first order: 



Ml = \ {gl + gl) 



2„,2 



= gj,L^' + -g'v 



(23) 



{gl + gl + g'Yv' - Kgl + gl - g')g%L^'x^ + le^^i^xM . (24) 



(25) 



All parameters in Eq. (16)-(20) as well as in the following mass matrices are understood as 



running parameters at a given renormalization scale Q. Note that the vevs Vd and f„ are 
obtained from the running mass Mz{Q) of the Z boson, which is related to the pole mass 
Mz through 



M|(Q) 



gt + g2 



{vl + vl) = Ml + Re{Ulz{Ml)}. 



(26) 



Here, is the transverse self-energy of the Z. See for more details also Ref. |71] . 



E. Neutralinos 



In the neutralino sector, the gauge kinetic effects lead to a mixing between the usual 
MSSM neutralinos with the additional states. Both sectors would decouple were these to be 
neglected. The mass matrix reads in the basis ( A^, H^, H^, \j^,,fj,f] 



TJly-O 
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9bl'^v 9bl^v 
-fx' 
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(27) 



It is well known that for real parameters such a matrix can be diagonalized by an orthogonal 
mixing matrix such that N*M^ N"^ is diagonal. For complex parameters one has to 
diagonalize Mf{Mf)^. 



F. Charged sleptons and sneutrinos 

We focus first on the sneutrino sector as it shows two distinct features compared to the 
MSSM. First, it gets enlarged by the superpartners of the right-handed neutrinos. Second, 
even more drastically, splittings between the real and imaginary parts of the sneutrinos occur 
resulting in 12 states: 6 scalar sneutrinos and 6 pseudoscalar ones [721 [73]. The origin of 
this splitting is the Y^^ OifjOj in the superpotential, Eq. ([T]), which is a AL = 2 operator 
after the breaking of U{1)b-l- Therefore, we define 

= ^ i^l + ^01) , ^1?, = ^ {^R + • (28) 

The 6x6 mass matrices of the CP-even (m?s) and CP-odd (^?p) sneutrinos can be written 
in the basis (ctl^o'r), respectively, {(Pl,4>r) as 

. S 1'"" V ™|,=»f'""'"- |. (29) 

^RL m^R / V m^RL ^RR 



While mj^]^ = = ttill hold^ the entries involving "right-handed" sneutrinos differ by 
a few signs. It is possible to express them in a compact form by 



rriLL = ml + "^YX + ^ {i9i + qI + f + 99bM - vl) + 2((7L + 99bl){'"1 ' ^-)) 1 



(30) 



m^l = (vuT* - VdfiY*^ ± VuV^Y:,Y* , (31) 



+ I {^91l{^1 - vD + -ggsL{< - vl)) 1 . (32) 

The upper signs correspond to the scalar and the lower ones to the pseudoscalar matrices 
and we have assumed CP conservation. In the case of complex trilinear couplings or /i terms, 
a mixing between the scalar and pseudoscalar particles occurs, resulting in 12 mixed states 
and consequently in a 12 x 12 mass matrix. In particular, the term ~ Vf^Y^fi'* is potentially 
large and induces a large mass splitting between the scalar and pseudoscalar states. Also 
the corresponding soft SUSY-breaking term ~ VrjT^ can lead to a sizable mass splitting in 
the case of large \A()\. 

The differences in the charged slepton sector compared to the MSSM are additional D 
terms as well as a modification of the usual D term. The mass matrix reads in the basis 
(cl, cr) as 

^ rriLL (vdTe - VulJ,*Ye 

rriLL = ml + ^Y^Y^ + ^ {{gl - gl + f + ggnM ' vl) + '^{99 bl + 91l){vI - 1 , 

(34) 

rriRR = ml + ^Y^Y} + ^ {{2gl + 2f + gg^,^) {vl -vl)-2 (2gg^^ + gl^) (t;^ - t;?)) 1 . 

(35) 

For the first two generations one can neglect the left-right mixing and thus the mass eigen- 
states correspond essentially to the electroweak flavor eigenstates. In the following we will 
call the partners of the left-handend (right-handed) leptons L sleptons {R sleptons). 



2 



ml 



(33) 



^ We have neglected the sphtting induced by the left-handed neutrinos as these are suppressed by powers 
of the light neutrino mass over the sneutrino mass. 
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G. High-scale boundary conditions 



We will consider in the following a scenario motivated by minimal supergravity, e.g. we 
assume a GUT unification of all soft SUSY-breaking scalar masses as well as a unification 
of all gaugino mass parameters 

2-1 2 2 2 2 2 /o'7\ 

Ml/2 =Mi = M2 = Ms = M^, . (38) 

Similarly, for the trilinear soft SUSY-breaking coupling the following mSUGRA conditions 
are assumed: 

Ti = AqYi, i = ejdjUjXjU . (39) 

Furthermore, we assume that there are no off-diagonal gauge couplings or gaugino mass 
parameters present at the GUT scale 

gsY =gYB = , (40) 
Mbb' =0 . (41) 

This choice is motivated by the possibility that the two Abelian groups are a remnant of a 
larger product group that gets broken at the GUT scale as stated in the Introduction. In 
that case gyy and gsB correspond to the physical couplings gi and g^i^, which we assume 
to unify with g2: 

gCUT ^ gCUT ^ gCUT ^^^^ 

where we have taken into account the GUT normalization discussed in Sec. Ill B[ 

In addition, we take the mass of the Z' and tan /3' as inputs and use the following set of 
free parameters 

mo, Ml/2, tan/3, tan/3', sign(yu), sign(/i'), Mz', Y^axidY^. (43) 

Yu is constrained by neutrino data and must therefore be very small compared to the other 
couplings, e.g. they are of the order of the electron Yukawa coupling. Therefore, they can 
be safely neglected in the following. Y^ can always be taken diagonal and thus effectively 
we have nine free parameters and two signs. 
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III. NUMERICAL RESULTS 



A. Implementation in SARAH and SPheno 

All analytic expressions for masses, vertices, RGEs, as well as one-loop corrections to 
the masses and tadpoles were calculated using the SARAH package [7iH76] . The RGEs are 
included at the two-loop level in the most general form respecting the complete flavor struc- 
ture using the formulas of Ref. [77] augmented by gauge kinetic mixing effects as discussed 
in Ref. [M]- The RGEs and the loop corrections to all masses as well as to the tadpoles are 
derived in the DR scheme and the Feynman-'t Hooft gauge. 

The numerical evaluation of the model is very similar to that of the default implemen- 
tation of the MSSM in SPheno [TSJ [TH]: as the starting point, the SM gauge and Yukawa 
couplings are determined using one-loop relations, as given in ref. [71], that are extended 
to our model. The vacuum expectation values Vd and Vu are calculated with respect to the 
given value of tan/3 at M^, while and Vf^ are derived from the input values of Mz' and 
tan /3' at the SUSY scale. 

The RGEs for the gauge and Yukawa couplings are evaluated up to the SUSY scale, where 
the input values of and are set. Afterwards, a further evaluation of the RGEs up to the 
GUT scale takes place. After setting the boundary conditions, all parameters are evaluated 
back to the SUSY scale. There, the one-loop-corrected SUSY masses are calculated using 
on-shell external momenta. These steps are iterated until the relative change of all masses 
between two iterations is below 10~^. 



We have used WHIZARD [80] to evaluate the bounds on the Z' discussed in Sec. Ill C as well 



as for the signals for slepton production in Sec. |III D[ For this purpose, the SUSY Toolbox 
[81] has been used to implement the model in WHIZARD based on the corresponding model 
files written by SARAH and to perform the parameter scans with SSP. 



B. Parameter studies 



In all numerical evaluations we have used the following SM input: Gp = 1.6639 • 
10-5 GeV-^ mz = 91.187 GeV, = 1.7771 GeV, rrit = 172.9 GeV, mb{mb) = 4.2 GeV 
and asimz) = 0.119. The latter two are MS values that are converted to the DR scheme. 
Moreover, we have fixed the neutrino Yukawa couplings such that the light neutrino masses 
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TABLE II. Parameters of the study points and selected masses. 
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FIG. 1. The mass spectra of BLV (left) and BLVI (right). 



can be explained. As this requires the maximum of \Xv,ii\ to be at most 10~^, these cou- 
phngs do not play any role in the considerations here. Because of their smallness, one can 
automatically satisfy the constraints from the nonobservation of rare decays such as /i — )■ 67 
over all the parameter space. This is in contrast to the usual MSSM augmented by the 
various seesaw mechanisms, as discussed in e.g. [821 183] and references therein. 

As a starting point for our numerical investigations, we have taken the two points defined 
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FIG. 2. Allowed regions in the mo - M1/2 plane for Mz' = 2.5 TeV, tan/3' = 1.15 (left) and in 
the Mz' — tan/3' plane for mo = 1 TeV and M1/2 = 1.5 TeV (right). In both plots we have fixed 
in addition tan/3 = 20 and Aq = —1.5 TeV. The grey dot indicates the selected benchmark point 
(BLV). 
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FIG. 3. Allowed regions in the mo — M1/2 plane for AIz' = 2 TeV, tan/3' = 1.07 (left) and in the 
Mz' — tan/3' plane for mo = 0.6 TeV and M1/2 = 0.6 TeV (right). In both plots we have fixed in 
addition tan/3 = 10 and Ao = 0. The grey dot indicates the selected benchmark point (BLVI). 
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in Table |TT| In analogy to Ref. [61], we label them BLV and BLVI (see Fig. Q. BLVI 
has a spectrum close to the existing LHC exclusion bounds. Note, however, that the high- 
scale input corresponds to different SUSY particle masses in comparison to the CMSSM 
because the running of the parameters is not the same and also the mass matrices differ. 
In contrast, BLV has quite a heavy spectrum and can only be discovered at ^/s = 14 TeV. 
As discussed above, in this model the extra gauge group implies additional D terms to the 
sfermion masses. The requirement that the mass-squared parameters for all the sfermions 
are positive restricts the allowed range for tan/3', as can be seen in Figs. [2]and|3} This is 
a consequence of the large mass for the Z'. Clearly this restriction is more severe for large 
Mz' and less severe for large values of mo and M1/2. 

For completeness, we note that these points would give too large a relic density, which 
however is not an insurmountable problem and can easily be fixed without changing the 
collider phenomenology. In case of BLV one would need to invoke nonuniversal boundary 
conditions for the bileptons to change its mass. In this way one obtains an efficient annihila- 
tion of the lightest neutralino via a bilepton resonance as discussed in [81] but having at the 
same time only tiny effects on the various decay branching ratios. In case of the BLVI the 
addition of a nonthermally produced gravitino with a mass of 10 GeV gives the correct relic 
density yielding a lifetime for the neutralino of about 10~^ seconds. This is sufficiently long- 
lived to appear as a stable particle at the LHC and at the same time sufficiently short-lived 
so that there are no problems with big bang nucleosynthesis. 



C. Z' phenomenology 

The mass of additional vector bosons as well as their mixing with the SM Z boson, which 
implies, for example, a deviation of the fermion couplings to the Z boson compared to SM 
expectations, is severely constrained by precision measurements from the LEP experiments 
[85H87] . The bounds are on both the mass of the Z' and the mixing with the standard model 



Z boson, where the latter is constrained by | sin(9vK') < 0.0002|. Using Eq. (22) together 



with Eq. (25) as well as the values of the running gauge couplings, a limit on the Z' mass 
of about 1.2 TeV is obtained. Taking in addition the bounds obtained from U, T and S 
parameters into account [HH] one gets Mz' / {Qf~'"gBL) > 6.7 TeV that for g^^ ~ 0.55 would 
imply Mz' > 1.84 TeV. However, this coupling has to be replaced by the effective coupling 
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FIG. 4. Branching ratios of the Z' at the considered parameter points BLV (above) and BLVI 
(below) as a function of the off-diagonal coupling parameter g. 



that gets modified due to gauge kinetic mixing: see Appendix k ' Therefore, the above 
formula reads in our model Mz' / {Q^^^^Qbl + d)) > 6.7 TeV, and g^^^ ~ 0.55, (? ~ —0.11 
imply Mz' > 1.47 TeV. 

The Z' dominantly decays into SM fermions as can be seen in Fig. |4] where we show the 
branching ratios as a function of g. We have fixed all other parameters as given for the 
^ Strictly speaking, this effect modifies the couplings to left- and right-handed fermions differently. Taking 

this into account would require a re-evaluation of the complete analysis which is beyond the scope of this 

paper. To be conservative we have taken the larger of the two couplings. 
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FIG. 5. LHC production cross sections of the Z' at the considered parameter points BLV (left) and 
BLVI (right) and three different center-of-mass energies as a function of the off-diagonal coupling 
parameter g. 



Parameter point 


g = - 


-0.11 


9 = 


= 


BLV 


1770 


GeV 


1965 


GeV 


BLVI 


1730 


GeV 


1900 


GeV 



TABLE III. Current bounds on Mz' in the supersymmetric B — L model derived from 5 fb ^ of 
ATLAS data for both benchmark points using different assumptions as discussed in the text. 



two study points in Table |TT} For the study points themselves we find g ~ —0.11, which 
has to be compared with g^^^ ^ 0.55. As they are of the same order of magnitude, one can 
easily understand the strong dependence of the various branching ratios after inspecting the 
couplings given in Appendix |A 1 The Z' can also decay into supersymmetric particles [53] 
with branching ratios of up to 0(10%) in our model. We find that, in particular, decays into 
charged sleptons can have a sizable branching ratio, as can also be seen in Fig. |4} Besides 
the decays, the cross sections also depend on gauge kinetic mixing as demonstrated in Fig. [5] 
where we show the Z' cross section for ^/s = 7, 8, and 14 TeV as a function of g. For the 
PDFs we have used the set CTEQ6L1 |89] . 

Recently ATLAS and CMS |44t US] have updated the results for Z' searches. As they 
do not give direct bounds for our model, we have calculated the corresponding signal cross 
section. If one took only final states containing SM particles into account and neglected 
gauge kinetic mixing, one would find a lower bound of 1970 GeV < Mz'- Taking into 
account gauge kinetic mixing we find for g ^ —0.1 a bound of 1790 GeV. Moreover, for 
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FIG. 6. Current limits on Mz' for combined lepton production for the study points BLV (left) 
and BLVI (right). The red dashed curve shows the recent experimental ATLAS limits [H]. The 
black and grey bands are the dilepton production cross sections at the Z' peak for the case of [/(I) 
mixing (black) and without (grey). The grey shaded area shows the mass range forbidden by LEP 
II, while the black dotted line shows the LEP limits without taking into account gauge kinetic 
mixing. 



both our study points, final states containing supersymmetric particles are also present as 
discussed above. This leads to an increase of the width and thus to a reduction of the 
signal cross sections. Therefore the bounds are less severe, as has also been discussed in 
the context of related models [56l [57] . In contrast to previous studies, we encounter here 
a case where neglecting gauge kinetic effects would lead to a significantly incorrect bound. 
In fact the bound obtained at LHC differs by about 200 GeV depending on whether kinetic 



mixing is correctly taken into account or not, as can be seen from Table III This is further 
exemplified in Fig. [6] where we display the cross section into leptons (summed over electrons 
and muons) as a function of Mz'- The grey area is excluded by precision LEP data whereas 
the red dashed line gives the bound on the signal cross section as obtained by the ATLAS 
Collaboration [H]. Note that the effect on the LEP limits is even stronger, as can be seen 
at the dotted line in Fig. [6j The CMS bounds are similar and thus lead to nearly the same 
limits. Clearly it would be desirable to have a combined analysis by both collaborations. 

In Fig. [rjwe give the cross sections for pp Z' at LHC with ^/s =14 TeV as a function 
of Mz' for the two benchmark points. To demonstrate the importance of gauge kinetic 
mixing, we display the values taking it into account (solid lines) and those where it is 
neglected (dashed lines). Comparing both data points, one sees that the dependence on 
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FIG. 7. Cross sections for pp — )• Z' at LHC with ^/s =14 TeV as a function of Mz' for the 
benchmark points BLV (left) and BLVI (right). The upper two hues show the total cross section 
whereas the lower two show the cross section in SUSY particles via a Z' . The solid lines give 
the cross section taking into account gauge kinetic mixing whereas the dashed ones give the cross 
section if gauge kinetic mixing were neglected. 

the gauge kinetic mixing on the SUSY final states depends on the underlying soft SUSY- 
breaking parameters. For BLVI, the dominant SUSY channels always involve the sleptons 
and therefore kinetic mixing reduces the cross section. In contrast, for BLV and small values 
of Mz'i final states including two neutralinos with large bileptino contents are important. 
Therefore, the cross section is larger with kinetic mixing than without. However, the masses 
of these neutralinos rapidly increase with Mz' and slepton production is the dominant SUSY 
channel for Mz' > 2.3 TeV. 

D. Slepton production via Z' as discovery channel 

A heavy Z' allows the production of electroweak SUSY particles with masses of several 
hundred GeV implying that this might be an important discovery channel as has also been 
discussed in related models In contrast to the previous studies, gauge kinetic mixing 

is important as we have seen above. In the following we will first take benchmark point BLVI 
to discuss some basic features. Much as in the MSSM, the R sleptons decay almost always 
into a lepton and the lightest neutralino, whereas the L sleptons decay dominantly into the 
lighter chargino and a neutrino, or the second lightest neutralino and a lepton, in roughly 
the ratio two to one. The lighter chargino and the second lightest neutralino, which have a 
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FIG. 8. Cross sections for the dominant final states resulting from the SUSY particles that are 
produced via a Z' for benchmark point BLVI. The symbols correspond to x = X?) 1 = (anti)lepton 
(e,/x), b = (anti)bottom quark, v = neutrino, j = jet resulting from a quark of the first two 
generations. 




FIG. 9. Feynman diagram for the production of II* via an s-channel Z' exchange and the subsequent 
decay into l^l^XiXi- 



large wino fraction, decay dominantly into the lightest neutralino and a vector boson. The 
neutralino decay can also result in a final-state Higgs boson. Therefore one has additional 
leptons and jets stemming from the decays of the vector boson and Higgs bosons. 

In Fig. [8] we display the most important final states resulting from the cascade decays 
of supersymmetric particles originating from the Z'. We have fixed the soft SUSY-breaking 
parameters as in BLVI. Nevertheless, the masses change as a function of Mz' as the corre- 
sponding bilepton vevs enter the mass matrices, see e.g. Sees. HE and |II F The dominant 



final state contains two leptons and two lightest supersymmetric particles stemming from 
slepton production as displayed in Fig. |9| For these points, sleptons are hardly ever produced 
in the cascade decays of squarks and gluinos, in contrast to the neutralinos and charginos. 
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FIG. 10. Histograms of the fJ-^n^ + If production with the apphed cuts M^^ > 200 GeV, pt[^) > 
200 GeV and Mj- > 800 GeV as a function of the missing energy, muon pair invariant mass and 
the transverse cluster mass. 
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For this reason the Z' decays are a potential discovery channel for sleptons. We have 
performed a basic Monte Carlo study using WHIZARD jHO] to generate the signal in the 
case of smuon pair production and simulated the background. For the background, we 
considered diboson production, triple vector boson production and tt production as well as 
neutralino and chargino production including both direct production via Drell-Yan processes 
and cascade decays from squarks and gluinos. However, the contributions from the latter 
are rather small. We have applied the following cuts to suppress the background 

• Invariant mass of the muon pair: M^^ > 200 GeV. 

• Missing transverse momentum: "Priff) > 200 GeV. 

• A cut on the transverse cluster mass 



Mt=^ + + - (pT(/i+/i-) +Pt(^))' (44) 

with "Pt being the 2D vector of the transverse momentum and pt(/^^/^ ) = Pr(yU^) + 
PT{fJ'~)- We required Mt > 800 GeV. 

For the suppression of ti and squark/gluino cascade decays, we set a cut on the trans- 
verse momentum of the hardest jet: prjet < 40 GeV. 



In Fig. 10 we display the resulting distributions for an integrated luminosity of 100 fb ^ at 



y/s = 14 TeV. The resulting significance, which is calculated as 

(45) 



A'signal 



is 7.5 a, which is sufficient to claim discovery for this example. As can be seen from 



Fig. 11, with a luminosity of about 45 fb~^ one crosses the 5a level. These numbers have 
been obtained using tree-level calculations and it is known that higher-order corrections are 
important in this context, see e.g. [90] and references therein. In case of nonsupersymmetric 
models one obtains K-factors of about 1.2-1.4 depending on the details of the models pUll^ . 
However, as it is not obvious how corrections due to supersymmetric particles will change 
this or how they affect the background reactions, we will stick to tree-level calculations here. 

The results depend so far mainly on the following quantities: Mz', mr, m^o and 
and on the nature of the neutralinos and charginos. We assume in the following that the 
lightest two neutralinos and the lighter chargino are mainly the MSSM gauginos as in our 
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FIG. 11. Evolution of the significance level with growing integrated luminosity. The discussed 
significance at 100 fb^^ fixes the curve. The borders for 3 and Scr are shown as dotted and dashed 
lines. 

study points. However, we will depart to some extent from the GUT assumptions: we will 
fix the masses of squarks and gluinos to the values they take in benchmark point BLVI in 
Table |Tl] but vary the slepton masses freely. As we are interested in relatively light sleptons 
with masses down to 200 GeV we fix m^o = 140 GeV and m^o = m-+ = 2m^o. We have 
performed a scan over slepton masses and Mz' fixing the ratio of the masses for R sleptons to 
L sleptons to 1.2, 1, and 1/1.2. The different ratios can in principle be obtained by varying 
tan as can also be seen from the formulas in Sec. II Fp We have also investigated the 
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effect of tightening the previous cuts to Pri^) > 250 GeV and M^^ > 300 GeV. In Fig. 
we show our results for the significance assuming an integrated luminosity of 100 fb~^. We 
find a significant dependence on the ratio of the smuon masses which is due to the fact that 
we consider here only smuon decays into fiXi- This is also the reason for the "promontory" 
for slepton masses up to 300 GeV as there both L and R smuons can decay only into this 
final state. The regions should even become somewhat larger if one includes also the smuon 
decays into fj,X2 ^Xi ■ However, one has then to consider the background of squark and 
gluino cascade decays that depend on the details of the parameter point under study. 



In Fig. 13 we display the same for the case of runj^ = 1.2 m/i^ but taking a luminosity 



^ Note, however, that tan /?' < 1 can only be obtained with non-universal bilepton mass parameters at the 
GUT scale. 
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FIG. 12. Significance level of the smuon discovery as a function of the respective lightest fl mass 
and Mz' for C = 100 fb^^ and m^o = 140 GeV. Left column: looser cuts [pt{^) > 200 GeV, 
M^f, > 200 GeV, Mt > 800 GeV]. Right column: tighter cuts [pt{^) > 250 GeV, M^^ > 300 GeV]. 
The smuon mass relations are (first row) m/i^ = 1.2 m^^, (second row) = m^^ and (last row) 
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FIG. 13. Significance level of the smuon discovery for C = 300 fb ^, rn^o = 140 GeV m^^ = 
1.2 rrijx^. Left: looser cuts [pt{^) > 200 GeV, M^^ > 200 GeV, Mt > 800 GeV]. Right: tighter 
cuts \pt{I^) > 250 GeV, M^^ > 300 GeV]. The mass of the lightest supersymmetric particle is 
140 GeV. 
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FIG. 14. Significance level of the smuon production for C = 100 fb ^, m^o = 280 GeV = 
1.2 m^^. Left: looser cuts [pt(^) > 200 GeV, M^^ > 200 GeV, Mr > 800 GeV]. Right: tighter 



cuts [pt{$) > 250 GeV, M^^ > 300 GeV] 



of 300 fb . As expected, the reach increases for both and Mz'- For the other two 



ratios of slepton masses we find the same behavior. In Fig. 14 we show the changes for 
increased neutralino and chargino masses m^o = 280 GeV and X2 = Xi = 475 GeV, which 
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are the values obtained in the constrained model at BLVI. As expected, the 5a significance 
is restricted to larger values of the smuon mass as the leptons are softer than in the previous 
example. 

IV. CONCLUSION 

We have studied a supersymmetric model where the gauge group is extended by an 
additional U{1)b-l factor resulting in a Z' with a mass in the TeV range. Such models can 
emerge as effective models from heterotic string models. An important feature of this class of 
models is gauge kinetic mixing. Here we focused on its impact for the phenomenology of the 
extra vector boson. We have first shown that bounds on its mass due to collider searches get 
significantly reduced once the gauge kinetic mixing is taken into account: LEP bounds by 
about 400 GeV and LHC bounds by about 200 GeV. Moreover, this imphes that now LHC 
bounds are more important than the ones originating from LEP data. Second we have shown 
that these bounds get further reduced if the Z' can decay into supersymmetric particles. In 
our model the most important decays are into sleptons, neutralinos and charginos. 

Moreover, we have discussed the reach of LHC with ^/s = 14 TeV for the discovery 
of smuons, which is an important possibility should squark or gluino cascades in sufficient 
quantities be inaccessible. For an integrated luminosity of 100 fb~^ (300 fb~^) we have 
found that sleptons with masses of up to 800 GeV (900 GeV) can be discovered this way, 
provided the Z' is lighter than about 2.8 TeV (3.1 TeV). This result depends only mildly on 
the nature of the smuons provided that their decays into muons are not suppressed. Apart 
from detector effects, similar results hold for selectrons. Staus, on the other hand, require 
a detailed study and we expect a reduced reach at the LHC due to the hadronic decays of 
the resulting r leptons. 
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Appendix A: Z' couplings 

In this section we collect the formulas for the couplings of Z' to the fermions and scalars 
in this model. 

1. Couplings to fermions 

The couplings given below follow from the terms in the Lagrangian 

]il^{ci^UiPL + c«/,,,Pi?)/,Z; . (Al) 

• Charged leptons: Z' — Cj — 

CLe,ij =\^ij ( (s'l sin - 9i cos ©ly) sin ©V + + 9bl) cos ©V) , (A2) 
CRe.ij = ^(^jj (^9\ sin sin ©V + (^9 + 9bl) cos ©V) ■ (A3) 

• Neutrinos: Z' — v^ — i/j 



CLu,ij =^ ( ( (f 1 sin + 92 cos Qw) sin ©V + + 9bl) cos ©V) ^ UjfU^a 

a=l 

3 

- 9bl cos ©V J] , (A4) 

a=l 

1 ^ 

Ci?,.,ij = ~ 2 ( ( (^^ ^'^ + + S-bl) cos ©V) UlfUj^ 

3 

-5BLCos©V5^^^t^]^+„). (A5) 

[/^ is the unitary 6x6 matrix that diagonalizes the neutrino mass matrix. 
Up-type quarks: Z' — Uia — Uj^ 

CLu,ij = - ^SaisSij (^(^-3g2 COS Qw + 9i sin ©u^) sin ©V + (§ + 9bl) cos ©V) , 

(A6) 

CRu,ij = - ^Sa^Sij (igi sin Qw sin ©V + (4=9 + 9bl) cos ©V) ■ (A7) 
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• Down-type quarks: Z' — dia — djp 

CLd,ij = - ^SapSij (3g2 cos + gi sin Qw) sin eV + (g + 9bl) cos ©V) , (A8) 
CRd,ij = ^SapSij (2gi sin Qw sin Q'w + (^^9 - 9bl) cos 6V) • (A9) 

• Neutralinos: Z' — x° — x° 

^Lx°,ij =^ (^i3 ( (91 sin ©ly + 92 cos 6^/) sin GV + 9 cos GV) A^is 

— A^*4 ^g^i sin ©vi^ sin ©V + 92 cos ©^ sin ©V + 9 cos ©V^ A^i4 

+ 2(?^^ cos G V [n*,N,, - N*,N„'^ ) , (AlO) 
CRxO,ij = - ^ (n^s ( (s'l sin ©1^ + g2 cos ©ly^ sin ©V + 9 cos ©V) -^^-3 

— iV*4 ^5(1 sin Qw sin G V + (72 cos Qw sin G V + 9 cos G -^^4 

+ '^9bl cos© [n*,Nj, - N:,Njr)) . (All) 

A^fci is the unitary 7x7 matrix that diagonahzes the neutrahno mass matrix. 

• Charginos: Z' - xj - xj 

CLx±,ij = - ^ (^92U*i cos Qw sin Q'wUn 

— U*2 (^gi sin Qw - g2 cos Qw^ sin ©V + 9 cos ©V^ Ui2^ , (A12) 
CRx^,ij = - ^ (^^92Vil COS ©H/ sin Q'wVji 

— V^2 ( (51 si^ — 92 cos ©H/^ sin ©V + 9 cos ©V^ Vj2^ ■ (A13) 

C/fe; and Vki are the unitary 2x2 matrices needed to diagonahze chargino mass matrix. 
2. Couplings to scalars 

The couphngs given below follow from the terms in the Lagrangian 

^s,i^s^~s]{p%-p':^|Z'^, (A14) 

where p^. and Ps*. are the four-momenta of the scalars. In the following, Z^^ denote the 
matrices needed to diagonalize the respective underlying mass matrix of the particles p. 
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Charged sleptons: Z' — Ci — e*j 

1 ^ 

Ce,y = 2 ( ( ^™ ~ + {9 + 9bl) cos OV) ^fa*^fa 

0=1 

3 

^2(71 sin Bh^ sin 9 V + {^g + <7bl) cos GV) Y1 ^^^ta^fs+a) ■ (A15) 



0=1 



Sneutrinos: Z' — Df — D^* 



-i^.y ( ~ ( (f 1 ®^ + f2 cos Qw) sin OV + + Qbl) cos OV) ^1 Zffzff 

a=l 

3 

- COS eV E ^£a^S;a) ■ (A16) 



Up-type squarks: Z' — — u* 



3 

f f — .Sr7o nns -I- rji sin (-)ti/^ sin f)\ir A- fn-\-n^r \ ms G'u/^ 



(^qu,ij = - (^(^(^- 3g2 COS + fifi sin sin GV + (5 + S'bl) cos GV) X] 



a=l 

3 

+ (4^1 Sin Qw sin G V + (4^ + 9bl) cos G V) ) ■ (A17) 

a=l 

Down- type squarks: Z' — di^ — d^p 

1 ^ 
-5c«/3 ( ( (35(2 COS Giy gi sin G^y) sin GV + (g + gBh) cos GV) X] Z^a* ^fa 



+ 



0=1 



- 2^1 sin Qw sin GV + ( - 2^ + ^^l) cos G V) ^f^aZf^^^ ■ (A18) 



a=l 



Charged Higgs: - //r - //+ 

CH±,ij ^\^i3 ( (fl'i sin - g2 COS GvK^ sin GV + g cos GV) ■ (A19) 

CP-odd and CP-even Higgs: Z' - - hj 

CAh,ij =^ ( - ( (s"! + 5*2 COS Qw^ sin GV + ^ cos GV) Z^Zf^ 

+ ^ (^gi sin Gvi/ + (72 cos Qw^ sin G V + g cos GV j ^^2^^^ 
- 2gs^ cos GV (^i^^i^ - Z^,Zf,) ) . (A20) 
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3. Coupling to vector bosons 



The only three- vector-boson vertex containing a Z' is the couphng Z'^ — — . It is 
parametrized as follows: 

cv.v, {g,, ( - pf" + p5) +9,.{- P^" + P^") +9.,{- P^" + pf^))z'^VfV[ , (A21) 
with 

cww = 92 cos Qw sin 6V . (A22) 

4. Coupling to one vector boson and one scalar 

The vertices are parametrized as follows: 

cvs,iSi9.^.Z'^V'' . (A23) 

• Z and Higgs: Z'^ — Z^, — hi 
2 



czKi = ^ ( ~ (fi'ifi' ®'vK sin 011/ + 92 cos 0^ cos 9 V sin 9V 



— (71(7 sin 9vy sin 9 



/2 

W 



COS 9V sin 9^ - (7^ j sin 9V - sin 
-I- g2 cos 9iy ^(71 sin 9vi/ sin 29V + ^7 cos 9'^ — (7 sin 9'^ j j 
— "^M {9x9 cos 9'^ sin 9vy + 9\ cos 9^ cos 9 V sin 9 V 
+ cos 9V {gX sin 9^ — (f^ sin 9V — fl'i^ sin 9vy sin 9'^ 
-I- g2 cos 9ty ^(^1 sin 9iy sin 29V + ^7 cos 9'^ — (7 sin 9'^ j j 
+ 2(7^^ sin 29V [vr^Zl + ) . (A24) 
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